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Global stability of traveling wavefronts in a periodic spatial-temporal environment in 
n-dimension (n > 1) is studied. The wavefront is proved to be exponentially stable in the 



. form of 0(e for some /i > 0, when the wave speed is greater than the critical one, and 

algebraically stable in the form of 0(t~"/^) in the critical case. A new and easy to follow 



method is developed. These results are then extended to the case of time-periodic media. 
Finally, we illustrate how the stability result can be directly used to obtain the uniqueness 
of the wavefront with a given speed. 



Keywords: Traveling wave solutions, stability, uniqueness, periodic media 

> 

^ ■ AMS 2000 MSG: 35K57, 92D25 

m ■ 
in 
cn 



1 Introduction 

Since the pioneering works of Fisher [TU] and KPP [T7] , wave propagations for reaction-diffusion 
models in fluid mechanics, physical, chemical and biological systems have been investigated in- 
tensively. For example, such wave behaviors with front propagation into unstable states can be 
frequently seen in fluid dynamical experiments on Taylor-CouettefH [19] and Rayleigh-Benard 
system [8], and in chemical wave experiment [281 131j as well as in population dynamics, com- 
bustion and biological invasions in a homogeneous environment [21 [3] . The theory of spreading 
speeds and monostable traveling waves has been developed to monotone semiflows in such a 
way that it can be applied to various evolution systems admitting the comparison principle 
(see |33[ [20l [TS] ) . For the stability of monostable waves including the minimal wave fronts in 
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some regular and time-delayed reaction-diffusion equations in spatially homogeneous habitat, 
we refer to [321 [T3l [22l [23| |25| [24] and the references therein. The stability for bistable waves 
to the regular or time-delayed reaction-diffusion equations has also been intensively studied by 
many people, for example, see [9l [211 HSl [ISl [301 132] and the references therein. 

Recently, the study of traveling waves and spreading speeds in heterogeneous media or 
periodically fragmented environments has been drawing more attention, see, e.g., the survey 
paper |30J. In a series of works, Berestyski et al. [H [H [7j and Hamel [12j established the 
existence of traveling fronts of the following equation 

ut-V ■ {A{x)Vu) = f{x,u), xGM", (1.1) 

where the diffusion matrix A{x) and the reaction term / both are periodic in x = (xi, X2, a;„) 
in the sense that there exist constants Lj > 0, 1 < "i < n, such that 

A{x + LiGi) = A{x), f{x + Liei,u) = f{x,u), VI < i < n, 

with ej being the i-th vector in the standard base of M" (i.e., its i-th component is 1 and all 
other components are 0). The matrix A{x) is also assumed to be symmetric and lies in the 
space C^''^{W^). The reaction term f{x,u) is a continuous function of class C^''^ with respect 
to X locally uniformly in u S M. The derivative is assumed to exist and be continuous in 
G M uniformly for x in M". For instance, / can be taken as f{x,u) = u{fi{x) — u), with the 
positive periodic function fi{x) usually being the carrying capacity of the periodic environment. 
Furthermore, for the matrix A, it is usually assumed to be bounded, and uniformly elliptic in 
the sense that 

^^A^MJ > «0|C|', "0 > (1.2) 

l<i,j<n 

for any vector ^ € M". 

In the special case where the period = 0, namely, the matrix A and the reaction term / 
are independent of x, model (jl.ip is a generalization of the simple Fisher-KPP model 

{ut — Au = fiu) 
^ ' (1.3) 
/(0) = /(l)=0, /(n) >0, Vne (0,1). 

It is well known that for any c > 2, system (II. 3p with f{u) = u{l — u) has a planar traveling 
wave solution u{t,x) = U{x ■ e-\- ct) with |e| = 1, where — e is the direction of propagation, 
while it has no such traveling wave for any c < 2. 

The existence of traveling fronts to (|l.ip has been studied recently in |11[T2]- For spreading 
speeds and traveling waves of other types of biological evolution systems in a periodic habitat, 
we further refer to [341 \TE\ [35] and the references therein. However, as mentioned in the survey 
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paper \30\ page 185], the stability analysis of traveling wavefronts in periodic media is a very 
challenging research topic. In particular, the author of |30j emphasized that the dynamics of 
the slow waves moving at the minimal speeds requires "a more delicate argument and awaits 
further investigation" . 

In this paper, we will concentrate on the stability and uniqueness of the wavefront of (jl.ip . 
an extension of the classical KPP-Fisher model in heterogeneous media. Motivated by the 
work of [26j on the classical Fisher-KPP equation and our recent study |24] on nonlocal time- 
delayed reaction-diffusion equations, we further investigate the stability of wavefronts for the 
periodic spatial-temporal reaction-diffusion equations by the Green function method. A new 
and easy to follow approach in terms of a wave transformation is developed. The wavefront is 
proved to be exponentially stable in the case where the wave-speed is greater than the minimal 
one (also called the critical speed), and algebraically stable in the case of the minimal speed. 
Furthermore, we extend our results to the case of time-periodic media. In the final concluding 
remark, we illustrate how the stability result can be directly used to obtain the uniqueness of 
the wavefront with a given speed. 

After we finished our first draft of this paper, we happened to obtain an important preprint 
|14j from Drs. Roques and Hamel's personal websites, where they investigated the stability and 
uniqueness of pulsating wavefronts for (jl.ip as well. It is seen that both papers are focusing on 
the open problem of the stability of wavefronts. However, there exist some major differences. 
Firstly, our developed method is quite different from that in |14j . Secondly, we prove that 
all non-critical pulsating wavefronts are exponentially stable in the form of 0(e"^*) for some 
^ > 0, and the critical pulsating traveling wave is algebraically stable in the form of 0{t^'"'^^), 
while there are no convergence rates in the global stability results of [14J. Thirdly, we assume 
that the product of the initial perturbation and the weight function is in L^{R^), while it is 
assumed in |14] that initial functions behave like a front. 

2 Exponential and algebraical stability 

In this section, we are going to investigate the global stability of planar traveling waves to 
the spatial periodic equation (jl.ip . Let w{x,x • e + ct) be a given planar traveling wavefront 
with wave speed c for some direction e ([e| = 1, — e is the direction of propagation) such that 
w{x + Liei,^) = w{x,^), Vx G M", ^ = x- e + ct£R, l<i<n. Without loss of generality, 
we can always assume that e = ei = (1,0, ...,0) after a rotation of coordinates. Note that the 
existence and the asymptotic behavior of w{x,S^), ^ = xi + ct, can be obtained from [H 112) . 

For equation (jl.ip with f{x,0) = to have wavefronts connecting two steady states, the 
following sub-linearity of the function / is also assumed in [U [12] : 
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s — )• ^ — is decreasing in s > for any x G M", (2-1) 

and 

f{x, s) < 0, s > M for some positive constant M. (2-2) 
We have the following simple observation. 

Lemma 2.1. // condition 12. 1\) holds, then we have 

fu{x,u)<Ux,Q), VxGR", n>0, (2.3) 

where fu{x,u) is the partial derivative of f with respect to the variable u. 

Proof. Define 

g{x,s) := Zi^, Vx €M", s > 0. 

Then ()2.ip implies that the function g is decreasing in s > 0, and hence, ^ < 0. Since 
f{x,s) = sg{x,s), we have 

df . . dg , . , . df(x,0) 
— = g{x, s) + s—< g{x, s) < g{x, 0) := • 

This completes the proof. □ 

The non- negative steady-states of p. II) . involving in the boundary condition of wavefront 
at iboo, are the C^'"(]R")-solutions p{x) defined by 

J V ■ {A{x)Vp{x)) + f{x,p{x)) = 0, 

\ p{x) > 0, p{x + LiGi) = p{x), Ml < i < n. 

Obviously, the constant is a trivial steady-state and its local stability is associated with an 
eigenvalue problem. Let /xi be the principal eigenvalue of the operator Lq defined by 

Lo<i)= V ■{A{x)V<i>) + fu{x,Q)(t>, 

that is, 

r V-(^(x)V0) + /„(x,O)0(x) =/ii0(x), 

{ (2.4) 

[ (t){x) > 0, (pix + LiGi) = </>(x), VI < i < n. 

If III > 0, then we say the zero solution is linearly unstable. In this case, there exists a unique 
positive and L-periodic steady state solution p{x) such that every positive solution u(t, x) of 
(|1.1|) with L-periodic initial function converges to p{x) uniformly in x G M" as t — > oo, see |5j. 

Now we can see that when fii > 0, there are two non-negative steady-states pi{x) = 
and P2{x) = p{x), which has the same mono-stability as in the Fisher-KPP equation. The 
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heteroclinic connection between pi{x) and P2{x) gives rise to a traveling wavefronts we are 
seeking. In other words, by a traveling wave solution, we mean a particular solution u{t, x) = 
'w{x,^),^ = xi + ct, which satisfies equation (jl.ip subject to the boundary condition 

lim w{x,s) = 0, lim w{x,oo) = p{x), uniformly in x G M". (2-5) 

S^ — OO S— )-CXD 

Substituting the relation u{t,x) = w{x,S^) with = xi + ct into (ll.ip . we obtain the wave 
profile equation 

r c^= V-{A{x)Vw{x,x,+ct)) + f{x,w), 

< (2.6) 

[ w{x, — oo) = 0, w{x, oo) = p{x), uniformly in x G M". 

Now we consider the case where > 0, that is, the steady state is unstable. The following 
result on the existence of spatially periodic traveling wavefronts is from ^tl2j. 
Theorem A ([11[T2]). Assume that /ii > and / satisfies (|2.ip and (|2.2p . Then for any c > c*, 
there exists a positive traveling wave u{t, x) = w{x, xi +ct) = wi^x,^) satisfying (jl.ip and (j2.5p . 
namely (12. 6p . where the minimal speed c* is determined by 

c* = min{c : 3 A > 0, such that ^c(A) = 0}. (2.7) 

Here ^c(A) is the principal eigenvalue of the operator 

- L^^x'ip = -V • (^(x)VV') - 2Ae- {A{x)Vxp) (2.8) 
- [AV • {A{x)e^) + A^e- {A{x)e^) - Ac + /„(x, 0)] V', 

that is, 

-Lc.aV' = IJ'cWi' with ^(x + LiBi) = ip{x), Vx G M", 1 < i < n. 

The derivation of c* in [U [12] is based on the technique developed by Aronson and Wein- 
berger [21 [3] in terms of a monotone-semiflow defined from equation (jl.ip . Here we can give 
an alternative, also easy to follow, explanation directly from the wave profile equation ()2.6p . 
This idea originates from study of the KPP-Fisher equation and now is called the technique of 
stability analysis of positive wavefronts near the unstable steady state in their phase plane. 

Due to biological or physical reason, we always expect a non-negative wavefront, and as 
such, near the far-field ^ = —oo we assume w ~ v{x)e^^ for some positive constant A and 
positive periodic function v{x), both of which are dependent on the wave speed c. Substituting 
this relation into ()2.6p and ignoring the higher terms, we then get 

-V-(^(x)Vu)-2Ae- {A{x)Vv)-[XV ■ (^(x)e^) + A^e- (A(x)e^) - Ac + /„(x,0)] v = 0, (2.9) 



which implies that the operator —L^^x defined in (12. Sp should have zero as its principal eigen- 
value. The minimal speed c is the smallest of such c that (j2.9p holds for A > 0. This agrees 
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with the definition in (|2.7p . Furthermore, in terms of A, the operator —Lc^\ is a polynomial 
of degree 2 and the result in [12] shows that the minimal speed c* is well-defined and has the 
following properties. 

Theorem B ([l2]). When c > c*, there exists two constants Ai(c),A2(c), < Ai(c) < A2(c) 
such that /.ic(Ai(c)) = /Xc(A2(c)) = 0, /Uc(A) > for any A € (Ai(c), A2(c)); wMe in the critical 
case c = c* , we have jic'W < 0, and ijLc*{^*) = when A* = Ai(c*) = A2(c*). 

To present the exponential stability of wavefronts, we also impose a natural condition on 
the positive steady-state p2 = p{x). We assume that it is linearly stable in the sense that the 
principal eigenvalue fii of the operator Li, defined by 



hold for some negative constant /xi and non-zero function 0. 

The subsequent result shows that the condition /ii < can be easily realized. 

Lemma 2.2. Assume that 

(C) fuix,u) < •^^'^'"^ , Vtt > 0, X € M", and there exists xo € M" such that fuixo,u) < 

Then any possible positive L-periodic steady state p{x) of is linearly stable. 

Proof. Let g{x,u) := f{x,u)/u, Vn > 0, and let fi{m{x)) be the principal eigenvalue of the 
operator Li with fu{x,p{x)) replaced by a continuous and L-periodic weight function m{x). 
Since p{x) = p{x + L) satisfies the elliptic equation and only the principal eigenvalue has 
positive eigenfunction, it is easy to see that fi{g{x,p{x)) = 0. Since fu{x,p{x)) < g{x,p{x)) 
and fuixo,p{xo)) < g{xo,p{xo)), the monotonicity of the principal eigenvalue with respect to 
the weight function implies that fi{fu{x,p{x))) < fj,{g{x,p{x))) = 0. Thus, p{x) is linearly 
stable. □ 

For any given wavefront w{x,S^), = xi + ct with speed c > c*, we define a weight function 



where is chosen in such a way that the wavefront w{x, ^), ^ > is very close to the positive 
steady-state p2 = p{x), that is, is chosen so that 



Li^ = V ■ {A{x)Vc^) + fu{x,p{x))<t), 



is negative, that is. 




(2.10) 




(2.11) 



w{x,i)-p{x)\ < e, Ve > Co 



(2.12) 
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uniformly in x E M" for any given sufficient small constant e > 0. 
Now we are in a position to prove the following stability result. 

Theorem 2.1. Assume that (j2.ip and (|2.'2|) hold, ni > 0, and fix < 0. Let u{t,x) be the 
solution U.l\) with the initial condition uq{x) satisfying < Uo{x) < p{x), and w{x,xi + ct) 
be a spatially periodic traveling wavefront with c> c*. If the initial data uq{x), having at most 
finite points of discontinuity, satisfies 

W{xi){uo{x) - w{x, xi)) G Li(M"), (2.13) 

then the following statements are valid: 

(1) For ^ = xi + ct < a holds 

\u{t,x) -w{x,xi + ct)\ < C(l + t)"te-'''=(^)*; (2.14) 

(2) For ^ = xi+ct> ^0, it holds 

\u{t,x)-w{x,xi + ct)\ < C(l + t)-te-"^'"^^^^(^)'^>*. (2.15) 

Here X, as mentioned in the weight function, is any number in the interval [Ai(c), A2(c)] . 

Proof. For the given initial function uo{x), we define 

I Uflix) = max{no(a;), w{x,xi)}, x € M", 
I Uq (x) = mm{uQ{x), w{x,xi)}, a; € M", 

where w{x,xi) is the initial function of the traveling wave solution w{x,xi + ct). 

Let U^{t,x) be the solutions of (jl.ip with the initial data U^{x), respectively, that is, 

^ = V- (A(x)VC/±) + /(x,?7±), , , 

< dt ^ ^ ^ ) J \ ^ )^ ^2.17) 

_ C/±(0,x) = U^{x), X G M". 
Since < < w{x,xi) < Uq < p{x), the comparison principle implies that 

< U'{t, x) < w{x, xi + ct) < U^{t, x) < p{x), 

and 

< U~{t,x) < u{t,x) < U+{t,x) <p{x). 

Obviously, we have 

U~{t, x) — w{x, xi + ct) < u{t, x) — w{x, xi + ct) < U~^{t, x) — w{x, Xl + ct), 
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or 



\u{t, x) — w{x, xi + ct)\ < max{|C/~''(t, x) — w{x, xi + ct)\, \U {t, x) — w{x, xi + ct)\}. 
As such, to prove (|2.14p and (j2.15p . we proceed to show 

\U^{t,x) -w{x,Xl + ct)\<C{l + t)-'^e-^'^^^^\ for ^ < -^o (2.18) 

and 

|[/±(t,x)-u;(x,xi + rf)l <C7(l + t)-ie-'°'"^^^=(^)'^>*, for ^ > Co- (2.19) 

In what fohows, we give a detailed proof of p.l8p and p.lOp for U^{t, x). The similar arguments 
work for U^{t,x), too, and hence, we omit them. 

Recall that both f7+(t,x) and w{x,xi + ct) satisfy the first equation of (|2.17p . Note that 
C = xi + ct. Thus we define 

V{t, x) = U^{t, x) - w{x, xi + ct). 
It is easy to see that V{t,x) is a nonnegative solution of the following equation 

Vt = V • (aw) +fix,w + V)- f{x, w) 
_y(0,x) = U^{x) -w{x,xi) > 0, 
with the initial data V{0,x) satisfying 

W{xi)V{0,x) = W{xi) {U^{x) - w{x,xi)) G L'^{W), (2.20) 

where W{xi) is the weight function defined in ()2.1ip . 

Since the non-linear function / is sub-linear and satisfies (j2.3p . it follows that 



f{x,w + V)-fix,w)<fu{x,0)V. (2.21) 

We then have an inequality 

'Vt<V-(AVv)+fu{x,0)V, 



(2.22) 

V{0,x) = U(^{x) -w{x,xi) > 0. 

Let v{x), dependent on c and A, be the eigenfunction of operator —Lc^x with the principal 
eigenvalue /Xc(A). From Theorem B, when c > c*, /^c(A) is real and v{x) > is well-defined. By 
a crucial wave transformation 

V{t, x) = 7;(x)e^(^i+^*-«»)y(t, x)e-'^=(^)*, A G [Ai(c), A2(c)], (2.23) 
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it then follows from (j2.22p that V{t, x) satisfies 

'Vt < V - (aV?) +git,x)-VV, 
V{0,x) = t;-i(x)e-^(^i-«o)F(0,x) > 0, 



(2.24) 



where 

g = 2A + Ae^ . (2.25) 
To get (12.24p . we have made use of the following formulas 

V{uv) = uVv + vVu, V • {vAVu) = vV ■ {AVu) + Vv ■ [AVu). 
It is easy to see that V in (I2.24p can be estimated by 

V{t,x)<l G{t,x -y)V{0,y)dy, yx£M.'\t>0, (2.26) 

where G{t,x — y) is the fundamental solution (i.e., the Green function) so that / G{t,x — 
y)uo{y)dy satisfies the following partial differential equation 

ut = V ■ {AVu) — g ■ S/u, n(0, x) = uo(x), 

with 

/ G{t,y)dy = l. 

Furthermore, since A is spatially periodic and uniformly elliptic, it follows from [11, (6.12) in 
page 24] that 

C i^-i/l^ 

\G{t,x-y)\<—e 4.1* (2.27) 

for some ai > 0, where n is the spatial dimension. By (j2.20p . we understand that V{0,x) G 
L^(M"). From ([226]), together with ([S^ZD, we can obtain 

Vit,x)<Gt-^ [ \V{0,y)\dy, Vx G M", t > 0. 

— n 

Since the solution V for t near zero has no singularity, the term t~2 can be replaced by 
(1 + t)~"/^. This technique is frequently used in the heat equation. Returning to (I2.23p . we 
then have 

V{t,x) = \U^{t,x) -w{x,xi +ct)\ 

= u(a;)e^(^'i+'=*-«o)F(t, x)e-^=W* 

< Ce^(^-i+^*-«o)(l + f)-fe-'^^(^)* (2.28) 
for ah X G M", t > 0. When ^ = xi + ct < ^o, we have e^^^"^") < 1. This proves (j2Tl]) . 
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Next for = xi + ct > (^q, we know that 

w{x,^) < w^x,^) + V< p{x). 

Since we have chosen so that is very close to p{x) uniformly in x € when ^ > 

(for precision, see (|2.12p ). as such for = xi + ct > by the regularity of /, we have 

fix, w + V)- fix, w) < ifuix,pix)) + e)V 

and 

Vt < V • iAVV) + ifuix,pix)) + e)V, 

ViO,x) = U^ix) - wix,xi) > 
for some small constant e satisfying < e < —fli, where fli is defined in (|2.10p with eigenfunc- 
tion (pix) > 0. Under the transformation 

Vit, x) = Vit, x)(j)ix)e^P^+'^\ (2.29) 

it then follows that for all ^ = xi + ct > ^(^j x) satisfies 

Vt<V ■iAVV) + 2(A^] - VV, ViO,x) = (j)-^ix)ViO,x), (2.30) 



with 

Vit,x) < y(t,x)(/.-i(x)e-(''i+^)* 
on the line xi + ct = ^o- Since on the boundary line ^ = xi + ct = (I2.28P yields 

Vit, x) < C(l + t)-te-^-W*, (2.31) 
which implies that on the line xi + ct = 

Vit,x) < C7(l + t)-te-'^^(^)*e^2*, 

where fj,2 = —fii — e > 0. Let 

/imin = mmj^c, ^-|- 
Thus, applying (j2.3ip to (|2.29p . we obtain that on the line xi + ct = 

Vit, x) < C(l + t)-te-''--*e^2t_ 

Obviously, C(l + t)~"2 e~'^™i"*e^^* is an upper solution of the equation 

ut = V ■ iAVu) + (2yl— -) • Vu, t>0 

4> 

for some constant C . By the parabolic comparison principle, it then follows that 
Vit,x) = 1/(t,x)(/)(x)e(^i+^)* < C7(l + t)-5e-''-'"*, Vx e M", t > 0. 
This completes the proof. □ 
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The following result is a straightforward consequence of Theorem 12.11 

Theorem 2.2. Assume that (j2.ip and (j2.2p hold, jii > 0, and fii < 0. Let u{t,x) be the 
solution lil.l]) with the initial condition uq{x), < ^0(2;) < p{x), and w{x, xi + ct) be a traveling 
wavefront with c > c* . If the initial data Uo{x), having at most finite points of discontinuity, 
satisfies 

W{xi){uo{x) - uj{x, xi)) e L\W), 

then the following statements are valid: 

(1) For any c > c* and A G (Ai(c), A2(c)), the wavefront w{x, x + ct) is exponentially stable 
in the sense that 

sup \u{t,x) —w{x,xi + ct)\ < Ce~^'^^, for some /X3 > 0; (2.32) 

(2) For c = c* , or c > c* with A = Ai(c), the wavefront is algebraically stable in the sense 
that 

sup \u{t,x) -w{x,xi +ct)\ < C(l + t)"t. (2.33) 

Proof For any given c > c* and A € (Ai(c), A2(c)), it follows from Theorem B that fJ,c{X) > 0. 
Let = min{ /ic(A), /i2}- Then we have (l2:32]l from (fmi) and (fZTSll . When c = c*, or c> c* 
with A = Ai(c), we then get /ic(A) = 0, and a combination of (I2.14p and (12.15P yields (I2.33p . □ 

3 An extension to the time-periodic case 

We have already studied the system which is periodic in space. In this section, we extend our 
results to the case when the system is also periodic in time. 
We consider the following parabolic equation 

ut-V ■ {A{t, x)\/u) = f{t, x,u), xe M", (3.1) 

where the matrix A and the nonlinear term / are periodic both in t and x, that is, 

A{t + T,x) = A{t, x), fit + T, X, u) = fit, X, u), 

and 

Ait, X + I/jej) = Ait, x), fit, X + Liei, u) = fit, x,u). Ml < i < n 

for some positive real number T > and n-dimensional vector L = (Li, L2, ...,Ln). 

As in [12, page 364-366], we assume that the uniformly elliptic symmetric matrix Ait,x) 
is of class Cl'"^'^'^'°' and the function fit,x,u) is continuous, of class (70'"/2;0,o respect to 
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(t, x) locally uniformly inii G R. We also assume that ^ exists and is continuous in u uniformly 
for X G M" and t G M . The existence of traveling waves w{t, x, xi + ct) was studied in |12l page 
364-366] and [27J with the sub-linearity assumption on the reaction term /. Assume that there 
are two periodic states pi{t,x) = and P2(t,x) > of (j3.ip . The system is called mono-stable 
if the state pi{t, x) is unstable and the other state P2it, x) is stable. In what follows, we assume 
that ^1 > and fli < 0, where is the principal eigenvalue of the linearized operator around 
0: 

ip{t, x) ^ ift — V ■ {A{t, x)Vip) — fu{t, X, 0)ip, ip{t, x) is periodic in t and x, 
and fli is the the principal eigenvalue of the linearized operator around p2{t,x): 

ip{t, x) — )■ (^t — V • {A(t, x)Vip) — fu(t, x,p2)^, (p{t, x) is periodic in t and x. 

Traveling wavefronts are solutions of (j3.ip in the form u{t, x) = w{t, x, e-x+ct) = w{t, x, xi+ct), 
where e = (1,0, ...,0), which is periodic in the first two variables t and x and satisfies 

lim 'w{t,x,^) = pi{t,x) = 0, lim wit, X, 6.) = P2itix) > 

uniformly for all t and x. The minimal speed c* is defined by 

c* = min{c : 3 A so that //c(A) = 0}, 

where n{X) is the principal eigenvalue of the operator 

v{t, x) ^vt + cXv-V ■ {A{t, x)S/v) - 2Ae- (AVv) - AV • (Ae) - A^e- (vie) - fu{t, x, 0)v, 

that is, there exists a periodic function v{t,x) > such that 

vt + cXv-V ■ {A(t, x)Vv) - 2Ae- (AVv) - AV • (Ae) - A^e- (^e) - f^it, x, 0)v = fic{X)v. 

It then follows that when c > c*, there exist < Ai(c) < A2(c) such that Hc{Xi) = ^Jic{X2) = 
and fJiciX) > for A G (Ai,A2)- In particular, when c = c*, we have Ai(c*) = A2(c*) and 
^c(A) < for all A. 

The weight function is defined as in (j2.1ip . The inequality (|2.22p now becomes 

Vt<V ■{AVV) + fi{t,x,Q)V, 
y(0, x) = U^{x) - w{Q, X, xi) > 0. 

Correspondingly, the transformation (|2.23p is replaced by 

Vit,C) = t;(t,x)e^(^-i+^*-«o)y(t,x)e-^^(^)*, A G [Ai(c), A2(c)]. 
All the remaining arguments in section 2 still work. Therefore, we have the following results. 
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Theorem 3.1. Let u{t, x) be a solution \3. 1\) with initial condition Uo{x), < uq{x) < P2it, x), 
and w{t,x,xi + ct) be a traveling wavefront with wave speed c > c* . If the initial data uq(x), 
having at most finite points of discontinuity, satisfies 

W{xi){uo{x) - w{0, X, xi)) G L\W), 

then the following statements are valid: 

(1) For ^ = xi + ct < it holds 

\u{t, x) - w{t, X, 01 < C{1 + t)-te~^<=W*; 

(2) For = xi + ct > ^,0, it holds 

\u{t,x) -w{t,x,0\ < C7(l + t)-?e-'"''^^^=W'^>*. 
Here A is any number in the interval [Ai(c), A2(c)] . 

Theorem 3.2. Let u{t, x) be a solution iS. 1\) with initial condition Uo{x), < ^0(2;) < ^2(^5 2;), 
and w{t,x,xi + ct) be a traveling wavefront with wave speed c > c* . If the initial data uo{x), 
having at most finite points of discontinuity, satisfies 

W{xi){uo{x) - w{0, X, xi)) G L\W), 

then the following statements are valid: 

(1) For any c > c* and A € (Ai(c), A2(c)), the wavefront w{t,x,xi + ct) is exponentially 
stable in the sense that 

sup |n(t, x) — w{t^ X, xi + ct) I < Ce~^^*, for some /U3 > 0; 

(2) For c = c* , or c > c* with A = Ai(c), the wavefront is algebraically stable in the sense 

that 

sup \u{t, x) — w{t, X, xi + ct)\ < C(l + t)~ 2". 

4 A concluding remark 

The stability analysis above, together with asymptotic behavior of spatially periodic wavefronts 
w{x, xi + ct) near the far field, can be directly used to obtain the uniqueness of the wave with 
a given speed. This then may provide a new and insightful approach that is different from the 
technique in |14j . 

To illustrate this idea, we first recall a result from [121 Theorem 1.3]. 



13 



Theorem C ([I2])- If O c* , then there exists a constant B > such that 



(4.1) 



where v{x) is the eigenfunction defined in ^2.9\) . If c = c* , then there exists a constant B > 
such that 

w{x,C) ^ B\C\e^^'^'''^^v{x), ^^-oo (4.2) 

where Ai(c*) = A2(c*) = A* and v{x) is the eigenfunction in i2. y\} with c = c* . 

Indeed, we can further claim that if c > c*, then the wave profile w{x,£,) has the following 
asymptotical behavior: 

w{x, C) = Se^i(^)^t;(x) + 0(e°^), B >0, a> Ai(c), ^ ^ -oo. (4.3) 

To see this, we may assume from (021) that w{x,^) = 5e^i W?v(x)(l+ w{x,^)). Substituting 
this into ()2.6p . we can get an equation for t(;(x,^). For this new equation, after a simple 
asymptotic analysis near ^ = — oo, we can find that w{x,(^) decays to zero exponentially as 
^ — )■ — oo. This confirms our claim. Now for any given c > c*, suppose that there are two 
spatially periodic wave profiles wi{x,^) and W2{x,^) of (jl.ip . A shift of distance s in W2 then 
enables W2{x, ^ + s) to share the same leading term in (j4.3p as that of wi. Let A in the weight 
function be Ai(c). Then when t = 0, we have e~'^^^'^^^^'^~^°\wi{x, xi) — W2{x, xi + s)) € L^(]R"). 
By Theorem 12. 2( it follows that 

lim sup \wi{x, xi + ct) — W2{x, xi + ct + s)\ = 0. 

Using a change of variable C = xi + ct, we then see that 

lim sup \wi({^ - ct,X2,--- ,Xn),0 - W2{{C - ct,X2,- ■ ■ ,Xn),^ + s)\ =0. (4.4) 

By the Li-periodicity of i(;i((xi, X2, • • • ,Xn),S,) and t(;2((xi, 0:2, ■ ■ ■ , x„), ^ + s) with respect to 
xi € M, it easily follows from (|4.4p that 

wi{x,C) = W2{x,S, + s), Vx = (xi,X2, • • • ,x„) e M", C G M. 

This shows that for any given c > c*, the wave profile w{x,^) is unique up to translation in ^. 
In the case where c = c*, the same idea applies, but more delicate details are needed. 

Acknowledgements. We are very grateful to two anonymous referees for careful reading and 
helpful suggestions which led to an improvement of our original manuscript. 
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